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INTRODUCTION 


Our  work  during  this  period  continued  to  have  as 
its  main  objective  the  study  of  the  thermodynamics  and  struc¬ 
ture  of  plasmas;  equilibrium  and  time  dependent.  We  used 
both  rigorous  analysis  and  approximation  schemes. 

In  addition  we  started  a  new  program  involving 
both  theoretical  analyses  and  computer  simulations  on  inter¬ 
acting  quantum  systems.  Our  first  publication  in  this  field 
will  appear  in  the  Rapid  Communication  section  of  Physical 
Review  in  February  1987.  We  have  also  made  progress  on  the 
statistical  mechanics  of  the  Zakharov  model  which  describes 
the  interaction  of  an  electromagnetic  wave  with  a  plasma.  The 
model  leads  to  a  nonlinear  set  of  partial  differential  equations 
which  resemble  in  some  ways  equations  arising  in  nonlinear 
field  theories  which  have  been  studied  by  statistical  mechanical 
methods.  We  have  obtained  some  very  exact  results  for  this 
model  in  one  dimension  and  are  now  carrying  out  some  computer 
simulations  on  it. 

The  general  nature  of  the  dynamics  describe  a  quantum 
system  in  interaction  with  a  heat  reservoir  is  also  being  investi 
gated  in  some  simple  model  system.  We  are  currently  writing  up 
a  report  on  this  work. 
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1.  Correlation  Functions  for  Nematic  Liquid  Crystals 
(Joel  L.  Lebowitz  and  John  W.  Perram) 

We  show  that  the  correlation  functions  of  a  system 
of  hard  ellipsoids  whose  axes  are  constrained  to  be  parallel  are 
simply  related  to  the  corresponding  functions  of  a  hard  sphere 
fluid  with  the  same  packing  fraction.  The  existence  of  simple, 
accurate  approximations  for  the  latter  system  make  the  oriented 
ellipsoid  fluid  a  suitable  "reference  system"  for  the  nematically 
ordered  phase  of  liquid  crystals  [1-3]  and  other  orientationallv 
ordered  molecular  fluids  such  as  polar  fluids  in  an  external  field. 

The  properties  of  such  fluids,  with  a  given  angular  distribution 
f(aj)  ,  may  then  be  obtained  via  "standard"  perturbation  theory. 

Methods  for  obtaining  f(oo)  near  the  perfectly  ordered  state  are 
discussed. 

(Jour.  Molec.  Phvs. ,  50,  6,  1207,  1983) 

2.  One-Dimensional  Models  of  Anisotropic  Fluids 
(Joel  L.  Lebowitz  and  Jerome  K.  Percus) 

Equilibrium  properties  of  fluids  of  approximately 
spherical  molecules  can  be  adequately  described  over  a  wide  range 
of  temperatures  and  densities  by  the  classical  statistical  mechanics 
of  point  particles  interacting  via  spherical  pair  potentials.  The 
situation  is  less  satisfactory  for  highly  nonspherical  molecular  fluids. 
It  is  the  aim  of  this  paper  to  seek  some  such  comparison  test  by 
considering  a  situation  in  which  the  fluid  is  confined  to  a  very 


narrow  channel. 
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3.  Sphericalization  of  Nonspherical  Interactions 
(Joel  L.  Lebowitz  and  Jerome  K.  Percus) 

A  density  and  temperature-independent  spherically  symmetric 
reference  potential  is  constructed  for  an  interacting  classical  fluid  of 
nonspherical  molecules.  It  annuls  the  first  order  correction  to  the  free 
energy  and,  in  special  cases,  the  second  as  well.  The  potential  is  a 
limiting  form  of  that  used  successfully  for  N9  by  Shaw  et  al. ,  and  re¬ 
produces  numerically  the  approximation  of  Rarbov  and  Gelbart  for 

homonuclear  dumbell  molecules. 

(J.  Chem.  Phys. ,  79(1) ,  1983) 

4.  Equivalent  Potentials  for  Equations  of  State  for  fluids  of 
Nonspherical  Molecules 

(G.  0.  Williams,  J.  L.  Lebowitz  and  J.  K.  Percus) 

We  analvze  the  extent  to  which  the  equation  of  state  and 
other  thermodynamic  properties  of  systems  of  hard  nonspherical  molecules 
can  be  obtained  from  a  density  independent  hard  sphere  reference  system. 
We  conclude  that  the  median  and  Barker-Fenderson  prescriptions  effective¬ 
ly  reproduce  all  data  now  available.  We  discuss  the  motivation  for 
these  two  formulations  in  detail. 

(J.  Chem.  Phys.  81(4) ,  1984) 


5.  Thermodynamics  of  Homonuclear  Diatomic  fluids  from  the  Angular  Median 
Potential 

(David  MacGowan ,  Eduardo  M.  Waisman,  Joel  L.  Lebowitz  and  Jerome  K. 

Percus) 

The  use  of  the  angular  median  potential  as  a  temperature- 
independent  spherical  reference  system  for  approximating  molecular  fluids 
is  tested  for  its  predictions  of  thermodynamics.  Calculations  have  been 
carried  out  for  a  wide  range  of  homonuclear  di atomics  with  continuous 
atom-atom  potentials  believed  to  be  representative  of  the  full  range  of 
simulation  data  available  for  such  systems.  The  results  for  the  pressure 
are  surprisingly  good  both  in  the  detonation  regime  and  around  the  triple 
point.  In  the  latter  case,  however,  the  internal  energies  for  highly 
elongated  molecules  with  attractive  potential  wells  are  considerably  too 
positive.  Comparison  with  other  perturbation  theories  indicates  that  the 
median  reference  system  gives  better  pressures  but  poorer  energies  than 
RAM,  and  that  in  many  cases,  especially  for  purely  repulsive  potentials, 
it  gives  results  of  comparable  accuracy  to  those  obtained  with  nonspheri¬ 
cal  reference  systems. 

(J.  Chem.  Phys.  80(6) ,  1984) 
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o .  The  Two-Dimensional  One-Component  Plasma  in  an  Inhomogeneous 
Background:  Exact  Results 

(A,  Alastuey  and  J,  L.  Lebowitz) 

We  study  the  general  inhomogeneous  two-dimensional  Jellium 
where  the  background  density  varies  in  one  space  direction  only.  At 
=  2,  explicit  functional  representations  of  the  one-  and  two-body 
densities  of  the  particles  are  derived  in  terms  of  the  electrostatic 
potential  created  by  the  background.  The  present  model  can  be  used 
for  describing  a  large  variety  of  charged  interfaces. 

(J.  Physique  4f5,  1859,  December  1984) 

7.  Low-Frequency  Electric  Microfield  Distributions  in  Plasmas 

(Carlos  A.  Iglesias,  Hugh  E.  DeWitt,  Joel  L.  Lebowitz, 

David  MacGowan  and  William  B.  Hubbard) 

We  evaluate  the  low- frequency  component  electric  micro- 
field  distribution  at  a  charged  point,  i.e.,  the  field  distribu¬ 
tion  due  to  ions,  possibly  of  different  species,  which  interact 
through  an  electron-screened  potential.  The  method  emploved 
is  an  adaptation  of  the  adjustable-parameter  exponential  (APEX) 
approximation  previously  developed  for  the  high-f requencv 
component  and  involves  a  noninteracting-quasiparticle  representa¬ 
tion  of  the  electron-screened  ions  designed  to  yield  the  correct 
second  moment  of  the  microfield  distribution.  The  APEX  results 
are  compared  to  Monte  Carlo  simulations,  and  we  find  good  agree¬ 
ment. 

(Phvs.  Rev.  A,  31.,  3,  March  1985) 
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8.  Long-Wavelength  Oscillations  in  an  Inhomogeneous  One-Component 
Plasma 

(Joel  L.  Lebowitz  and  Ph.  A.  Martin) 


The  perfect  screening  of  charge  fluctuations  in  an 

equilibrium  plasma  is  extended  to  the  time-displaced  structure 

function  of  a  general  inhomogeneous  one-component  plasma.  We  find 

that  the  long-wavelength  modes  oscillate  undamped  with  a  single 

frequency  u.u2  being  an  angular  average  of  squares  of  plasma 
2  2 

frequencies  or  ■  4i re  p/m  in  uniform  systems  with  density  p.  Our 
results  are  derived  rigorously  from  the  Bogoliubov-Born-Green-Kirkrwood- 
Yvon  hierarchy  under  some  reasonable  assumptions  on  the  spatial  decay 
of  correlations  and  contain  as  special  cases  previously  obtained  re¬ 
sults  of  this  kind. 

(Phys.  Rev.  Letts.,  54,  14,  8  April  1985) 

9.  The  Ideally  Polarizable  Interface:  Integral  Equations 
(M.  L.  Rosinberg,  L.  Blum  and  J.  L.  Lebowitz) 


The  integral  equations  used  in  the  microscopic  theory  of 
the  electric  double  layer  were  extended  to  the  case  of  an  impermeable 
interfacr  separating  two  conducting  media  (ionic  solutions  or  plasma). 
This  system  is  a  model  for  an  ideally  polarizable  interface.  Exact 
relations  are  given  for  the  contact  values  of  the  one  particle  density 
function,  and  also  for  the  pair  correlation  functions.  We  solve 
numerically  the  Poisson-Boltzmann  (PB),  the  hypernetted  chain  (HNC) , 
and  mean  spherical  (MSA)  approximations,  and  compare  the  results  to 
the  exact  solution  of  the  one  component  plasma  in  two  dimensions. 

(J.  Chem.  Phys.  83(2),  15  July  1985) 


10.  Time-Dependent  Correlations  in  an  Inhomogeneous  One-Component 
Plasma 

(B.  Jancovici,  J.  L.  Lebowitz  and  Ph.  A.  Martin) 

Sum  rules  describing  perfect  screening  at  equilibrium  in 

a  classical  plasma  were  extended  to  the  time  displaced  structure 

function  of  an  inhomogeneous  one-component  plasma.  We  found  that 

there  are  long-wavelength  modes  which  oscillate  undamped  with  a 

—  —2 

single  frequency  u  ui  being  an  angular  average  of  the  squared 
plasma  frequency  at  infinity.  Our  results  were  derived  heuristically , 
allowing  also  for  quantum  effects,  from  linear  response  theory,  and 
rigorously  from  the  classical  BBGKY  hierarchy  under  some  reasonable 
assumptions  on  the  spatial  decay  of  correlations.  Special  cases 
were  investigated,  in  particular  plasmas  bounded  by  walls  of  varied 
shapes. 

The  sum  rules  express  the  macroscopic  property  of 
screening  which  is  possessed  by  Coulomb  fluids.  They  can  be 
derived  by  using  linear  response  theory  and  assuming  that  an 
external  charge  introduced  in  the  fluid  induces  a  polarization 
charge  which  cancels  the  external  one.  More  rigorous  derivations 
using  the  BGY  hierarchy  can  be  given  under  suitable  clustering 
assumptions. 

In  the  present  work,  we  derive  a  dynamical  and  quantum- 
mechanical  generalization  of  sum  rules,  for  inhomogeneous  one- 
component  plasmas  (OOP).  The  restriction  to  OCP's  (i.e.  to 
systems  of  identical  particles  of  charge  e  and  mass  m  in  a 
fixed  background  of  opposite  charge)  has  to  be  made  because  we 


use  a  unique  property  of  these  systems:  the  long-wavelength  plasma 
oscillations  are  not  damped.  From  a  microscopic  point  of  view, 
this  is  related  to  the  validity  of  a  dipole  sura  rule.  We  can 
however  deal  with  a  large  class  of  inhomogeneous  OOP's:  we  rIIow 
the  background  density  pb(q)  *  a  function  of  the  position  q  ,  to 
have  a  value  at  infinity  which  may  depend  upon  the  direction  n 
in  which  q  recedes  to  infinity.  More  precisely,  lira  pb(r,fi)  = 
exists  (for  almost  every  fi) ,  with  q  ■  (r,fJ),  r  »  |q|,  n  = 
angles  of  q.  This  allows  in  particular  for  OOP  plasmas  bounded  by 
walls  of  varied  shapes. 

Defining  the  usual  dynamical  structure  factor  S(q1#  t  |  q) 
a£  the  correlation  function  between  the  charge  density  at  time  t 
and  point  q^^  and  the  charge  density  at  time  0  and  point  q,  we 
find  the  proper  dynamical  and  quantum-mechanical  generalization  which 


reduces  to  the  classical  limit,  and  to  the  quantum-mechanical 
static  limit;  the  appearing  in  these  equations  Is  an 


angular  average  of  the  plasma  frequency  (squared) 
infinity. 


P 


at 


Note  that  a  system  of  real  atoms  and  molecules  in 


which  the  nuclei  are  treated  as  fixed  is  included  In  our  scheme  - 
provided  that  they  are  confined  to  a  bounded  region  of  space. 

This  region  is  then  imagined  surrounded  by  a  more  smeared  out 
background  pb(r,0)  which  has  the  limit  pm(Q)  as  r  ■*  *  . 

(Jour.  Stat.  Phys.,  41,  December  1985) 
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11.  A  Solvable  Model  for  Localized  Adsorption  in  a  Coulomb  System 
(M.  L.  Rosinberg,  J.  L.  Lebowitz  and  L.  Blum) 

A  model  for  an  Interface  with  localized  adsorption  is 
presented,  in  which  the  surface  has  a  distribution  of  sticky 
adhesive  sites  in  contact  with  a  Coulomb  fluid.  Contrary  to  the 
current  literature  on  the  electrical  double  layer  the  surface 
charge  is  in  dynamic  equilibrium  with  the  bulk  fluid.  The  sum- 
rules  obeyed  by  the  one-and  two-body  correlation  functions  are 
investigated.  Explicit  results  are  obtained  for  a  solvable  model, 
the  two  dimensional  one-component  plasma  at  reduced  temperature  2. 

The  effect  of  the  granularity  of  the  adsorbed  charge  on  the  ad¬ 
sorption  isotherm  is  discussed. 

When  a  solid  material  such  as  an  electrode  plate  or  a  colloid 
particle  Is  in  contact  with -an  electrolyte  solution  or  a  plasma,  it 
can  acquire  an  electrical  charge  by  specific  adsorption  of  ions. 

This  adsorbed  surface  charge  modifies  the  distribution  of  the  ions 

and  the  molecules  in  the  vicinity  of  the  surface  and  induces  a  sur¬ 
face  potential.  In  the  case  of  chemiaorbed  systems  (or  strong 
phya 1 aorpt ion)  the  adsorbed  ions  are  localized  on  adsorption  sites 
with  a  distribution  which  reflects  in  some  way  the  structure  of  the 
adsorbent,  as  for  Instance  a  periodic  lattice.  It  la  not  an  easy 
task  to  include  this  surface  charge  which  is  in  dynamic  equilibrium 
*ith  the  electrolyte  in  a  statistical  mechanical  treatment  of  the 
Interface.  for  this  reason,  in  almost  all  works  devoted  to  the 
study  of  CoulcxaJb  systems  near  a  wall,  the  surface  charge  is  a  pre¬ 
scribed  'external*  charge,  uniformly  distributed  on  the  wall. 

In  particular,  the  adsorbed  charges  cannot  'react"  to  the  bulk 
system. 


-  9  - 

In  OUT  work  we  found  an  exact  aolution  to  a  two 
dimensional  nodal  of  the  true  adaorptlon  problem.  The  Coulomb 
fluid  la  the  classical  one-component  plaaroa  (OCP)  at  the  special 
temperature  T  for  which  the  coupling  constant  T  3  6e2  -  2, 

S  “  (KgTJ  \  and  the  specific  (short  range)  adsorption  potential 
ia  the  sticky  potential  of  Baxter  . 


The  one  component  plasma  is  the  simplest  model  of  a 
Coulomb  fluid  and  recently  various  physically  Interesting  in¬ 
homogeneous  situations  have  been  worked  out  exactly  in  two 


dimensions  for  r  ■*  2 


.  They  are  restricted,  however,  to 


the  case  where  the  background  density  or  the  external  potential 
varies  only  in  one  direction. 

We  shall  consider  here  the  following  three  situations  : 
i)  a  semi-infinite  plasma  in  the  vicinity  of  an  impenetrable  wall 
(l.e. ,  a  line)  with  continuous  adsorption  }  11)  an  infinite  plasma 
in  the  vicinity  of  discrete  adsorption  sites  regularly  spaced  on 
a  line  j  Hi)  a  semi-infinite  plasma  in  the  vicinity  of  an  impe- 
netrable  wall  with  discrete  adsorption. 

while  i)  is  a  simple  extension  of  the  results  obtained 
in  (3,4,5),  it)  and  ill)  are  real  two  dimensional  problems  and 
the  structure  of  the  fluid  (for  instance,  the  one-  and  two-body 
distribution  functions)  will  vary  both  in  the  directions  perpen¬ 
dicular  and  parallel  to  the  line  of  sites. 

The  method  that  we  use  for  solving  this  problem  involves 

Infinite  expansions  in  terms  of  the  n-body  correlation  functions 

of  the  unperturbed  fluid  (l.e.  without  adsorption  sites).  While 

this  solution  is  valid  for  any  fluid  subject  to  an  external  one- 

body  potential  of  Baxter's  type,  this  is  so,  however,  only  in  a 

formal  sense  since  the  n-body  correlation  functions  are  generally 

unknown  quantities.  The  two-dimensional  OCP  at  T  "  2  is  a  rare 

exception  :  the  gaussian  structure  of  the  correlation  functions 

(4,11)  makes  the  infinite  expansions  summable. 

(Jour.  Stat.  Phys. ,  44,  1986) 
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12.  Monte  Carlo  Simulation  of  Hard  Spheroids 

(J.  W.  Perram,  M.  S.  Wertheim,  J.  L.  Lebowitz  and  G.  0. 
Willi ams ) 

We  present  Monte  Carlo  simulations  of  the  equation 
of  state  and  radial  distribution  function  for  a  model  fluid 
composed  of  hard  spheroids. 

(C hem.  Phvs.  Letts.,  105 ,  3,  1984) 


13.  Exact  Results  for  the  Two-Dimensional  Two-Component  plasma 
(D.  Nicolaides) 

The  free  energy  and  correlation  functions  of  a  two- 

dimensional  two-component  plasma  are  explicitly  computed  at  the 

2 

special  temperature  T  =  e"/2k  . 

u  B 

(Phvs.  Letts.,  103A,  1,  1984) 

14.  Exact  Results  for  the  Two-Dimensional  One-Comnonent  Plasma 
(D.  Nicolaides) 

We  show  that  the  free  energy  and  correlation 

O 

functions  of  the  2D  OCP  at  the  special  temperature  TQ  =  e"/2kR 
can  be  explicitly  computed  by  solving  a  free-field  Dirac  equation. 

(Phvs.  Letts.,  103A,  5,  1984) 


-il¬ 


ls.  Fluids  with  Highly  Directional  Attractive  Forces.  I. 
Statistical  Thermodynamics 

(M.  S.  Wertheim) 


A  new  formulation  of  statistical  thermodynamics  is 
derived  for  classical  fluids  of  molecules  that  tend  to 
associate  into  dimers  and  possibly  higher  s-mers  due  to 
highly  directional  attraction.  A  breakup  of  the  pair 
potential  into  repulsive  and  highly  directionally 
attractive  parts  is  introduced  into  the  expansion  of 
the  logarithm  of  the  grand  partition  function  in 
fugacity  graphs.  The  bonding  by  the  directional  attraction 
is  used  to  classify  the  graphs  and  to  introduce  a  topo¬ 
logical  reduction  which  results  in  the  replacement  of  the 
fugacity  by  two  variables:  singlet  density  p  and  monomer 
density  .  Results  for  the  thermodynamic  functions  as 
functionals  of  p  and  p^  are  given  in  the  form  of  graph 
sums.  Pair  correlations  are  analyzed  in  terms  of  a  new 
matrix  analog  of  the  direct  correlation  function.  It  is 
shown  that  the  low  density  limit  is  treated  exactly, 
while  major  difficulties  arise  when  the  Mayer  expansion, 
which  employs  only  p,  is  used.  The  intricate  resummations 
required  for  the  Mayer  expansion  are  illustrated  for  the 
case  where  dimers  are  the  only  association  products. 

(Jour.  Stat.  'Phvs.,  35,  19,  1984) 


Fluids  with  Highly  Directional  Attractive  Forces.  II. 
Thermodynamic  Perturbation  Theory  and  Integral  Equations 

(M.  S.  Wertheim) 

The  formalism  of  statistical  thermodynamics  developed 
in  the  preceding  paper  is  used  as  a  basis  for  deriving 
tractable  approximations.  The  system  treated  is  one  where 
repulsion  and  highly  directional  attraction  due  to  a  single 
molecular  site  combine  to  allow  the  formation  of  dimers, 
but  no  higher  s-mers.  We  derive  thermodynamic  perturbation 
theory,  using  the  system  interacting  with  only  the  repulsive 
potential  as  a  reference  system.  Two  distinct  integral 
equations  for  the  pair  correlation  are  derived.  The  first 
one  treats  both  parts  of  the  interaction  approximately; 
the  other  one  employs  the  repulsive  reference  system  used 
in  perturbation  theory.  We  show  that  each  of  these  integral 
equations  permits  the  calculation  of  an  important  thermo¬ 
dynamic  function  directly  from  the  solution  at  a  single 
state  of  density  and  temperature.  In  the  first  case  this 
applies  to  a  pressure  consistent  with  the  compressibility 
relation,  in  the  second  to  the  excess  Helmholtz  free 
energy  over  the  reference  system. 

(Jour.  Stat.  Phys. ,  35,  35,  1984) 
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17.  Sum  Rules  for  Inhomogeneous  Coulomb  Fluids,  and  Ideal  Conductor 
Boundary  Conditions 

(B.  Jancovici) 

The  charge  correlation  function  of  an  inhomogeneous 
fluid  was  shown  to  obey  a  sum  rule  involving  electrical  multioole 
moments  of  arbitrary  order;  this  sum  rule  is  a  further  generalization 
of  the  Sti Hinge r- Lovett  second-moment  condition.  An  application 
is  given  for  a  fluid  bounded  by  an  ideal  conductor  wall;  it  is 
shown  how  the  chemical  potentials  control  both  the  bulk  and  surface 
properties . 

(Submitted  to  J.  Phys. ,  France) 


18.  The  Dipole  Phase  in  the  Two-Dimensional  Hierarchical  Coulomb 
Oas:  Analyticitv  and  Correlations  Decay 

(0.  Benfatto,  0.  Callavotti  and  F.  Nicolo) 

We  illustrate  the  mechanism  nroducing  the  dipole  phase 

in  a  two  dimensional  Coulomb  system  by  a  detailed  analysis  of  a 


hierarchical  model.  We  prove  the  analyticitv  of  the  pressure  and 

n  o 

of  the  correlations  for  a  =  Se^  >  8tt  (i.e.  right  above  the 


usually  conjectured  value  for  the  Kosterlitz-Thouless  phase  transi¬ 
tion).  We  find  also  a  power  law  decay  for  the  correlations  with 


exponent  ct^/2Tr  as  the  hierarchical  distance  goes  to  infinity. 
(Commun.  Math.  Phys.  106 ,  277,  3986) 


